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Abstract. We discuss the possibility of obtaining model-free bounds on volatil- 
ity derivatives, given present market data in the form of a calibrated local 
volatility model. A counter-example to a wide-spread conjecture is given. 
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1. Introduction 

"... it has been conjectured that the minimum possible value of an option on 
variance is the one generated from a local volatility model fitted to the volatility 
surface.'" ; Gatheral Gat06l page 155]. 

Leaving precise definitions to below, let us clarify that an option on variance 
refers to a derivative whose payoff is a convex function / of total realized vari- 
ance. Turning from convex to concave, this conjecture, if true, would also imply 
that that the maximum possible value of a volatility swap (f(x) = x 1 ' 2 ) is the one 
generated from a local volatility model fitted to the volatility surface. Given the 
well-documented model-risk in pricing volatility swaps, such bounds are of imme- 
diate practical interest. 

The mathematics of local volatility theory (a la Dupire, Derman, Kani, ...) is 
intimately related to the following 

Theorem 1 ( Gy686j). Assume dY t = fi (t, lj) dt+a (t, u>) dB t is a multi- dimensional 
Ito-process where B is a multi- dimensional Brownian motion, fj,, a are progressively 
measurable, bounded and aa T > e 2 I for some e > (<j t denotes the transpose of 
a). Then 

(1.1) dY t = Viae (t, f t ) dt + a loc (t, F t ) dB t , Y = Y , 

Vioc(t,y) = E [/j, (t, u) | Y t = y] , 

a loc (t,y) = E[a(t,co)a T {t,uj)\Yt = y^ , 

(where the power \ denotes the positive square root of a positive definite matrix) 

has a weak solution Y t such that Y t l == Y t for all fixed t. 

We will apply theorem [1] only in the simple one dimensional (resp. two dimen- 
sional in section 4 below) setting where it is well known that the solution to (jl.lj) 
is unique (cf. |Kry67| or [SV061 Capter 7]). 

A generic stochastic volatility model (already written under the appropriate 
equivalent martingale measure and with suitable choice of numeraire) is of the 
form dS = SadB where a — a (t, u) is the (progressively measurable) instantaneous 
volatility process. (It will suffice for our application to assume a to be bounded 
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from above and below by positive constants.) Arguing on log-price X — log S 
rather than S, 

(1.2) dX t = a (t, uj) dB t - (a 2 (t, u) /2) dt, 

a classical application of theorem Q] yields the following Markovian projection re- 
suliQ: the (weak) solution to 

(1.3) dX t = a loc (t, a) dB t - (ofoc (t, X t ) /2j dt, X = X , 

(a loc (t,x) = E[a 2 (t,Lo)\X t = x] i ) 

has the one-dimensional marginals of the original process X t . EquivalentljQ, the 
process S — expX, 

dS t = aioc{t, S t ) S t dB t , 
known as (Dupire's) local volatility model, gives rise to identical prices of all Eu- 
ropean call options C(T,K)E It easily follows that af oc (t,S) is given by Dupire's 
formula 

(1.4) <(r .S)| s=K = 2^| s . 

Volatility derivatives are options on realized variance; that is, the payoff is given 
by some function / of realized variance. The latter is given by 

V T := (log 5) T = (X) T = f a 2 {t,uj)dt; 



o 



in the model dS — a (t, lo) S dB and by 

V T := (log S) T = (X) T = [ a 2 loc {t, X t ) dt 

Jo 

in the corresponding local volatility model. 

Common choices of / are / (ar) = x, the variance swap, f (x) = x 1 / 2 , the volatility 
swap, or simply f (x) = (x — K) + , a call-option on realized variance. See [FG05 
for instance. As is well-known, see e.g. |Gat06| . the pricing of a variance swap, 
assuming continuous dynamics of S such as those specified above, is model free in 
the sense that it can be priced in terms of a log-contract; that is, a European option 
with payoff log St- In particular, it follows that 



E 



E[V t 



Of course this can also be seen from (|1.3I) , after exchanging E and integration over 
[0,T]. Passing from Vr to f{Vr) for general / this is not true, and the resulting 
differences are known in the industry as convexity adjustment. We can now 
formalize the conjecture given in the first lines of the introductiorQ. 

Conjecture 1. For any convex f one has ElfCVr)] < S[/(Vr)l- 

Our contribution is twofold: first we discuss a simple (toy) example which pro- 
vides a counterexample to the above conjecture; secondly we refine our example 
using a 2-dimensional Markovian projection (which may be interesting in its own 
right) and thus construct a perfectly sensible Markovian stochastic volatility model 
in which the conjectured result fails. All this narrows the class of possible dynamics 



^Let us quickly remark that Markovian projection techniques have led recently to a number of 
new applications (see IPit06l . for instance). 

2 The abuse of notation, by writing both o"; oc (f, X) and oi oc (t, S), will not cause confusion. 

■^We emphasize that C(T,K) denotes the price at time t = of European call with maturity 
T and strike K. 

%t is tacitly assumed that / (Vr) , /(Vr) & re integrable. 
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for S for which the conjecture can hold true and so should be a useful step towards 
positive answers. 

2. Idea and numerical evidence 

Example 2. Consider a Black-Scholes "mixing" model dS = SadB,So = 1 with 
time horizon T = 3 in which a 2 (t,u>) is given by a\{t) or cr^(t), 

[2 if te [0,1], (2 if te [o,i\, 

al(t):=i3 tftG]l,2], a 2 _{t):=\l tfte]l,2], 
(l tfte]2,3], [3 ift€]2,3], 

depending on a fair coin flip e — ±1 (independent of B). Obviously V = V3 = 
J Q a 1 dt = 6 in this example, hence E[(V — 6) + ] = (V — 6) + = 0. On the other 
hand, the local volatility is explicitly computable (cf. the following section) and one 
can see from simple Monte Carlo simulations that for V = V3 



E (V-6)- 



0.026 > 



thereby (numerically) contradicting conjecture^ with f (x) = (x — 6) + . 

Our analysis of this toy model is simple enough: in section [3] below we prove that 
P[V = 6] 7^ 1. Since E[V] = E[V] = 6 and (x — 6) + is strictly convex at x = 6, 
Jensen's inequality then tells us that E[(V - 6)+] > = E [(V - 6)+]. 

The reader may note that an even simpler construction would be possible, i.e. 
one could simply leave out the interval [0, 1] where cr+ and a^L coincide. We decided 
not to do so for two reasons. First, insisting on cn_(t) = o-_(t) for t £ [0, 1] leads to 
well behaved coefficients of the SDE describing the local volatility model. Second, 
we will use the present setup to obtain a complete model contradicting conjecture 
[T] at the end of the next section. 

3. Analysis of the toy example 

We recall that it suffices to show that V = f af oc (t, Xt) dt is not a.s. equal to 
V = 6. The distribution of X t is simply the mixture of two normal distributions. 
More explicitly, X t = I {e=+1} X t , + + J {e= _ 1} X t 

X t ,±= J^a ± {s)dB s -\ J^a 2 ± {s)ds ~ ^(|E±(t),E±(i)), 
where E±(t) := j£ cr±(s) ds. Thus o-f oc (t,x) = E[a 2 (t, u)\X t = x] is given bsfl 

g +W Pxr , r _ (^+S+(t)/2) 2 1 , °-(t) r _ (x+S-(t)/2) a -| 



(3.1) af oc (t,x) 



1 exD T _ 1 r _ ( a! +S-(t)/2)- 

^/^y ex PL 2S+(t) J + ^/sT^^L aE_(t) 



Since (Ti oc = o~i oc {s,x) is bounded, measurable in t and Lipschitz in x (uniformly 
w.r.t. t) and bounded away from zero it follows from [SV06L Theorem 5.1.1] that 
the SDE 

dX t = a loc (t, dB t - \a 2 loc (t, X^j dt 



^More general expression for local volatility are found in IBM06I Chapter 4] and |Lee0 1 , HL09 . 
Note the necessity to keep a 1 (. , u>) constant on some interval [0, e] , for otherwise the local volatility 
surface is not Lipschitz in x, uniformly as t — > 0. 
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has a unique strong solution (started from Xg — 0, say). Since ai oc is uniformly 
bounded away from it follows that the process (X t ) has full support, i.e. for every 
continuous ip : [0, 3] — > R, <p(0) = and every e > 

P[\\X t -ip(t)\\ oo . [0tT] <e]>0. 

Indeed, there are various ways to see this: one can apply Stroock-Varadhan's sup- 
port theorem, in the form of |Pin 95, Theorem 6.3] (several simplifications arise in 
the proof thanks to the one-dimensionality of the present problem); alternatively, 
one can employ localized lower heat kernel bounds (a la Fabes-Stroock [FS86j ) or 
exploit that the Ito-map is continuous here (thanks to Doss-Sussman, see for in- 
stance [RW001 page 180]) and deduce the support statement from the full support 
of 5. 




Figure 1. Time evolution of local variance af oc (t,x) in depen- 
dence of log-moneyness. The bright strip indicates a set of paths 
with realized variance strictly larger than 6. 

Figure 1 illustrates the dependence of of oc (i, x) on time t and log-moneyness x. 
To gain our end of proving that V (cj) = J Q 3 crf oc (t, X t ) dt is not constantly equal 6, 
we can determine a set of paths {X t {uj)) for which V is strictly larger than 6. In view 
of Figure 1 it is natural to consider paths which are large, i.e. X t (uj) 6 [8, 10], for 
t e]l, 2 — i] and small, i.e. |X t (w)| < 1, on the interval ]2, 3]. A short mathematica- 
calculation reveals that V (uj) > 6.65 > 6 for each such path and according to the 
full-support statement the set of all such paths has positive probability, hence V is 
indeed not deterministic. 

Using elementary analysis it is not difficult to turn numerical evidence into rig- 
orous mathematics. Making (|3.ip explicit yields that of oc (i, x) = 2 for t € [0, 1] and 
that 

(2x+2+3t) 2 (2x+2+tf_ 
3 „ 8(2+3t) i 1 ~ S(2+t) 
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(2x+5+t) 2 (2x+3+3t) 2 
8(5+t) i 3 ~ 8(3+3t) 



(3.3) a f oc{t + 2,x) = ^± t 



(2x+5+ty (2x+3+3ty 
8(5+t) i 1 c ~ 8(3+3t) 



/5+t J V3+3T 

i i ■ , W . t 1 i . \ t lilt! , . 



for t €]0, 1]. We fix e G]0, 1] and observe that it is simple to see that lim^^co of oc (i+ 
l,a;) = 3, uniformly w.r.t. t G [e, 1], and that lim^^o °f oc (^ + 1> X ) > 2, uniformly 
w.r.t. t e]0, 1]. It follows that there exists some S > such that 

°foc(* + 1j x ) > 3 — e for a; > t e [e, 1] and 
CT z 2 oc (t + l,a;) > 2-e for |x| < d,te]0,l]- 

Thus we obtain 

(3.4) V(u) = / crf oc (t, X t (u;)) dt > 1 • 2 + (1 - 2e) ■ (3 - e) + 1 • (2 - e) 



for every path X{ui) satisfying X t {ixi) > | for t e [1 + e, 2 — e] and |X t (a;)| < (5 for 
t € [2, 3]. This set of paths has positive probability and the quantity on the 

right side of (I3.4[) is strictly larger than 6 provided that e was chosen sufficiently 
small. Hence we find that V is not constantly equal to 6 as required. 

For what it's worth, the example can be modified such that volatility is adapted 
to the filtration of the driving Brownian motion. 

The trick is to choose a random sign e, P(e = +1) = P(e = —1) = | depending 
solely on the behavior of (B t )o<t<i and in such a way that S\ is independent of e. 
For instance, if we let m(s) be the unique number satisfying P{S\/2 > m(s)\S\ = 
s) = P{S\/2 < m(s)\S\ = s) = \, it is sensible to define i := +1 if Si/ 2 > m(Si) 
and e := — 1 otherwise. 

We then leave the stock price process unchanged on [0, 1], i.e. we define a 2 (t) = 
a 2 (t) = 2 and S t = S t for t G [0, 1]. On ]1, 2] resp. ]2, 3] we set & 2 (t) := 2 + e resp. 
<7 2 (i) :— 2 — e and define 5t, t G]l, 3] as the solution of the SDE 

(3.5) dS t = a{t)StdB u S l =S l . 

Here (|3.5H depends only on S\ and the process (B t — £>i)i<t<3; since both are 
independent of e, we obtain that (St)i<t<3 and (St)i<t<3 are equivalent in law. It 
follows that V = J Q 3 a 2 (t,(jj) dt = 6 and since 5*4 and St have the same law for each 
t € [0,3], they induce the same local volatility model and in particular the same 
(non deterministic) V. 

4. Counterexample for a Markovian stochastic volatility model 

Recall that X denotes the log-price process of a general stochastic volatility 
model; 

dX t = a (t, u) dB t - (a 2 (t, ui) /2) dt, 
where a = a (t, uj) is the (progressively measurable) instantaneous volatility pro- 
cess0. Recall also our standing assumption that a is bounded from above and below 
by positive constants. We would like to apply theoremQ]to the 2D diffusion (X, V) 
where dV = a 2 dt keeps track of the running realized variancfl We can only do so 
after elliptic regularization. That is, we consider 

dX t = a(t,u)dB t - (a 2 (f,u;)/2) dt, 

da% = a 2 (t, u) dt + e 1/2 dZ t 



^We note that we intend to insert later the volatility of examplefS] but so far our considerations 
hold in general. 

^In other words, Vt = Jg T o" 2 (t, oS) dt, if Vq = which we shall assume from here on. 
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where Z is a Brownian motion, independent of of the filtration generated by B and 
a. It follows that the following "double- local" volatility model 

d±t = a dloc (t, XI, of) dB t - (a 2 dloc (t, XI of) /2) dt, 

da-t = °l loc {t,Xl,~a\)dt + e x l 2 dZ t , 

(with a\ loc (t, x, a) = E[a 2 (t, co) \X t = x, a\ = a]) has the one-dimensional marginals 
of the original process (X t ,af). That is, for all fixed t and e, 

, - law , ~r law p 

X t = X\ and a t = a\. 

Let us also note that the law of a\ is the law of V t = a° convolved with a standard 
Gaussian of mean and variance e. The log-price processes X and X s induce 
the same local volatility surface. To this end, just observe that X t = Xf implies 
identical call option prices for all strikes and maturities and hence (by Dupire's 
formula) the same local volatility: 

af oc (t,x) = E[a 2 (t,u)\X t = x] = E[a 2 dloc (t, X^ a s t ) \X s t - x\. 

Since the law of a time inhomogeneous Markov process is fully specified by its 
generator, it follows that the law of the local volatility process associated to (X) 
has the same law as the local volatility process associated to (X s ). 

We apply this to the toy model discussed earlier. Recall that in this example, 
with T = 3 

V T = I <J 2 {t,uj)dt = 6 
Jo 

whereas realized variance under the corresponding local volatilty model, 

V T = j\f oc (t,X t )dt 

was seen to be random (but with mean Vt, thanks to the matching variance swap 
prices). As a particular consequence, using Jensen 




= (E J a 2 (t,uj)dt-6y = (V T ~6)+ =0. 

We claim that this persists when replacing the abstract stochastic volatility model 
(X) by (X £ ), the first component of a 2D Markov diffusion, for any e > 0. Indeed, 
thanks to the identical laws of the respective local volatility processes the left- 
hand side above does not change when replacing (X) by the local volatility process 
associated to (X e ). On the other hand 

e[ <j 2 dloc (t,Xl,~al)dt = E(h T -e l ' 2 Z T ) 
Jo 

= E (a £ T ) — E (a T ) 

= E(V T +e 1/2 Z T ) = V T - 

Thus, insisting again that the process X is (in law) the local volatility model asso- 
ciated to the double local volatility model (X £ ,a £ ) we see that 

c = E (J o <r 2 oc {t,X t )dt-Q) + > (eJ a 2 dloc {t,X £ t ,~a £ t )dt-Q) + = 0. 
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(Observe that c > is independent of e.) Using the Lipschitz property of the 
hockeystick function, again Gyongy and the fact that a € T is normally distributed 
with mean Vt and variance sT we can conclude that 



G 



E^J a 2 dloc (t,Xl~at)dt-^f = E^J aj loc (t, X%, a s t )dt + 

< E((a E T - 6)+ + \siz T \) 

= E(a £ T - 6)+ + E\e^Z T \ = 3y/eT/2ir. 

Now we choose e small enough such that 3yJeT/2Tr < c, whence the conjecture fails 
to hold true in the double local volatility model for e > small enough. 

5. Conclusion 

Summing up, the double-local volatility model constitutes an example of a con- 
tinuous 2D Markovian stochastic volatility model, where stochastic volatility is a 
function of both state variables, in which conjecture[T]fails, i.e. in which the minimal 
possible value of a call option is not generated by a local volatility model. 
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